All-Sky Convolution for Polarimetry Experiments 
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We discuss all-sky convolution of the instrument beam with the sky signal in polarimetry experi- 
ments, such as the Planck mission which will map the temperature anisotropy and polarization of 
the cosmic microwave background (CMB). To account properly for stray light (from e.g. the galaxy, 
sun, and planets) in the far side-lobes of such an experiment, it is necessary to perform the beam 
convolution over the full sky. We discuss this process in multipole space for an arbitrary beam 
response, fully including the effects of beam asymmetry and cross-polarization. The form of the 
convolution in multipole space is such that the Wandelt-Gorski fast technique for all-sky convolu- 
tion of scalar signals (e.g. temperature) can be applied with little modification. We further show 
that for the special case of a pure co-polarized, axisymmetric beam the effect of the convolution 
can be described by spin-weighted window functions. In the limits of a small angle beam and large 
Legendre multipoles, the spin-weight 2 window function for the linear polarization reduces to the 
usual scalar window function used in previous analyses of beam effects in CMB polarimetry experi- 
ments. While we focus on the example of polarimetry experiments in the context of CMB studies, 
we emphasise that the formalism we develop is applicable to anisotropic filtering of arbitrary tensor 
fields on the sphere. 
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I. INTRODUCTION 

Over the past decade a number of increasingly so- 
phisticated experiments have reported detections of the 
temperature anisotropy in the cosmic microwave back- 
ground (CMB). Following the detection of degree-scale 
anisotropy by the COBE satellite [1], ground-based and 
balloon-borne experiments have pushed back the limits 
on resolution and sensitivity to provide estimates of the 
anisotropy power spectrum up to Legendre multi- 
poles of I w 700 (see e.g. Ref. [2] for a review of the situ- 
ation pre-BOOMERANG [3,4] and MAXIMA-1 [5,6]). 
The anisotropy power spectrum encodes a wealth of cos- 
mological information (e.g. Ref. [7]) in a highly com- 
pressed form, making it a very convenient data product 
from which to determine cosmological parameters (see 
Ref. [8] and references therein). 

The combination of Thomson scattering and the non- 
zero temperature quadrupole, as the radiation begins to 
free stream through recombination, leads to the robust 
prediction that the CMB should be linearly polarized, 
with an r.m.s. level of a few percent of the temperature 
anisotropics [9-12]. Detection of polarization would pro- 
vide complementary information to that obtained from 
temperature measurements, e.g. the unique polarization 
signature of gravitational waves (and vector modes) pro- 
vides the best hope of detecting their presence at last 
scattering [13,14]. Currently, only upper limits exist 
on the degree of linear polarization (e.g. Refs. [15-17]; 



see also Ref. [18] and references therein for a recent re- 
view) but detections should be made with the MAP [19] 
and Planck [20] satellites, several experiments from the 
ground [21-24], and the flights of the MAXIPOL [25] and 
enhanced BOOMERANG [4] balloon experiments. 

In an ideal linear polarimetry experiment, a given de- 
tector is configured to respond only to a single compo- 
nent of the electric field of the incident radiation, along 
incident directions contained within a small solid angle 
Ail (the beam size). For such an ideal experiment, the 
detector measures the flux (/ -I- Q)Ail/2 where / (to- 
tal intensity) and Q are the Stokes parameters of the 
incident radiation along the beam direction, and the po- 
larization basis vectors have been chosen with the x di- 
rection aligned with the polarimeter. In practice, such 
ideals are never achieved, and for precision polarimetry 
experiments, e.g. Planck, it is essential to take full ac- 
count of several beam effects on the measured signal. The 
polarization on the sky must be convolved with the re- 
sponse pattern of the detector, which will not be perfectly 
axisymmetric. Furthermore, the system will generally 
have some cross-polar contamination (non-zero response 
to more than one polarization component). Often the 
beam response pattern has non-negligible far side-lobes 
which are highly polarized due to reflection and diffrac- 
tion effects in the instrument, making a full-sky convo- 
lution necessary if the effects of stray light from bright 
regions (such as the sun, moon and galaxy in CMB exper- 
iments) are to be properly accounted for. The implica- 
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tions of a subset of these effects for the analysis pipeline 
of total power (unpolarized) experiments have been stud- 
ied recently in Ref. [26]. 

In this paper we present a formalism that allows all 
such non-ideal effects to be taken account of exactly and 
efficiently in multipole space. The form of the convolu- 
tion for polarized data in multipole space is very similar 
to that for unpolarized radiation. This fact allows the 
recently suggested algorithm of Wandelt and Gorski [27] 
to be used to compute rapidly the detector output for an 
arbitrary pointing direction and orientation of the detec- 
tor. The formalism presented here should prove useful in 
simulation and modelling of precision polarimetry exper- 
iments, as well as in the actual analysis of experimental 
data. Within the multipole formalism, it is simple to in- 
voke approximations, such as axisymmetry of the beam, 
where appropriate, to reproduce the approximate results 
used in previous analyses of beam effects in CMB po- 
larimetry experiments [28,29]. 

The paper is arranged as follows. In Sec. II we discuss 
the derivation of the detector response from its far-field 
radiation pattern, and introduce the beam response ten- 
sor. The multipole expansion of the beam in spherical 
scalar and tensor harmonics is presented in Sec. III. The 
scanning of the detector on the sky is described by a 
(time-dependent) rotation of the beam from some stan- 
dard configuration. This rotation is most conveniently 
handled in multipole space, and is described in Sec. IV, 
where an efficient algorithm for performing the rotation 
and convolution is also described. Section V deals with 
the case of an axisymmetric beam. It is shown that for 
a certain geometry of the polarized beam response, it is 
possible to describe the convolution of the beam with the 
incident linear polarization in terms of spin-weight 2 win- 
dow functions, first introduced in Ref. [29]. We give the 
window function for Gaussian beams of arbitrary angular 
size, and show that in the small-angle limit we reproduce 
the results of Ng and Liu [29]. We summarise our dis- 
cussion in Sec. VI. An appendix provides some details 
on the rotation properties of the tensor harmonics which 
are needed for Sec. IV. 



II. FAR-FIELD RADIATION PATTERN AND 
THE BEAM RESPONSE 

It is convenient to characterise the response of the de- 
tector and feed system by the far-field radiation pattern 
it would emit if used as a transmitter, rather than a re- 
ceiver [30]. Assuming a (quasi-)monochromatic system, 
the electric field in the far field is of the form 

E oc -SR{f exp[z(fcr - ojt)]}, (1) 
r 

where £ is a complex, transverse vector function on the 
sphere. (We refer to geometric objects having no com- 
ponents outside the surface of the sphere as being trans- 
verse.) Here r is radial distance and u = ck is the (mean) 



angular frequency of the radiation, where k 1/r is the 
wavenumber. The most general detector and feed system 
will produce a partially polarized signal in transmission, 
so that £ will generally be a slowly varying function of 
time (compared to w). However, only the stationary sta- 
tistical properties of £ are important for determining the 
response of the system to incident radiation. With the 
system in some specified orientation, the power received 
dH^tot when illuminated by the sky along some direction e 
is proportional to the intensity in that polarization mode 
which is the time reverse of Eq. (1). It follows that 

dWtot oc {\£ ■ £\^)dn, (2) 

where £ is the complex representative (or analytic signal) 
of the incident electric field propagating along — e, dfl is 
the clement of solid angle, and angle brackets denote time 
averaging over the slow variations in £ and £. Writing 
the components of the fields £ and £ on the orthonormal 
basis vectors erg and cr^ of a spherical polar coordinate 
system as e.g. £g and the contribution to Wtot can be 
written as 

^^cx^{li + QQ + UU-VV), (3) 

where {I,Q,U,V} are the Stokes parameters of the in- 
coming radiation on the {crg,—cr^} basis: 

I={\£g\' + \£^\'), (4) 

Q={\£g\'-\£^\'), (5) 

U = -2^{£g£;), (6) 

V = 2^{i£g£;), (7) 

and {/, Q, U, V} are effective Stokes parameters for the 
beam: 

i={\£gf + \£^\'), (8) 
Q = {\£g\'-\£^\'), (9) 
U = -21k{£g£l), (10) 

V = 2^{i£g£;). (11) 

Note that the effective Stokes parameters for the beam 
are defined on the same basis as the incoming radiation, 
which is responsible for the minus sign in front of the VV 
term in Eq. 3. 

The intensity / and circular polarization V are in- 
variant under rotations of the polarization basis vec- 
tors, whilst Q and U transform like the components of a 
second-rank tensor. Introducing the linear polarization 
tensor for the incident radiation 

V^e) = ^[QicTg ®(Tg-(T^® (T^) 

-U{(Tg®(T4, + (T^®(TB)], (12) 

we can write the total power received in the basis- 
independent form 
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Wtot <x^J{li-VV + 2VabB'''>) dO. (13) 
Here, Bab is the (linear) beam response tensor: 

-U{cr0^cr^ + cr4,^ ere)] (14) 
= 5?[(£0r)rT, (15) 

where '^'^ denotes the transverse, trace-free part. Our 

main task now is to derive the dependence of the total 
power received on the pointing direction and orientation 
of the detector. 



A. Co- and cross-polarized basis 

The polar basis {erg, cr^} is fixed relative to the sky and 
is singular at the north and south poles. For describing 
the beam, it is standard practice to use an alternative 
basis which is fixed relative to the detector, and has its 
only singularity in the opposite direction to the main 
beam [31]. We define a set of Cartesian basis vectors 
{<T^, cr^, cr^} which are fixed relative to detector. It is 
convenient to take cr'^ to be along the (nominal) main 
beam, and cr'y along the polarization direction on axis. 
Using this Cartesian frame we derive a set of polar basis 
vectors {Tg,cr'^} on the sphere in the standard manner. 
The CO- and cross-polar basis vectors are then derived 
by parallel-transporting cr'y and cr^ respectively from the 
north pole along great circles through the poles: 

(Tco = sin (/)' ctq + cos (/)' cr^ (16) 
f cross = cos (j)' a'g — sin (j)' cr'^, (17) 

where 9' and (j)' arc spherical polar coordinates. A well- 
defined linearly polarized receiver has |f cross | \£co\ 
along the main beam, where e.g. ^co is the component 
of £ along ctco- Cross-polar contamination arises from 
a non-zero | Across |- To rotate from the co- and cross- 
polarized basis to the spherical polar basis we have to ro- 
tate through 71/2 — 4)' in a right-handed sense about the 
inward normal to the sphere. Transforming the Stokes 
parameters for the beam from the co- and cross-polar 
basis to the spherical polar basis, we have: 



I — (i^col + l^crossl )■, (18) 

Q=-(|5co|^-|^crossP)cOs2<^' 

+ 25J(fco4*ross)sin20', (19) 

i7=-(|£co|'-|^cross|')sin20' 

-2SR(fco4*ross)cos20', (20) 

y=-23?(z^co4*ross), (21) 



which reflect the spin-2 nature of the linear polarization. 
For simulation purposes, £ is usually determined with 



physical optics codes, the results being reported on the 
CO- and cross-polar basis. For real experiments the Stokes 
parameters for the beam must be calibrated using sources 
with known surface brightness and polarization. Note 
that we do not assume that the beam response is fully 
polarized, so the formalism developed here can also be 
applied to total power experiments (only I non-zero). In 
practice, the optics and feeds will introduce some beam 
polarization in the sidc-lobcs even for a nominal total 
power experiment. Although small, the role of such ef- 
fects in total power experiments could be quantified with 
our formalism. 



III. MULTIPOLE EXPANSIONS 

The dependence of the total power received on the di- 
rection and orientation of the telescope is most easily 
formulated in multipole space. We say that the detector 
is in its reference orientation when it is oriented so that 
the basis {cr'^,a'y,a'^} coincides with the {(Tx,(Ty,(Tz} 
basis, which is fixed relative to the sky. We describe the 
beam via a set of constant multipole coefficients which 
are extracted on the sky basis when the detector is in its 
reference orientation. To describe an arbitrary orienta- 
tion of the detector at some time along the scan, we can 
rotate the beam (which is most easily performed in mul- 
tipole space) to obtain the rotated beam response which 
is convolved with the sky, as in Eq. (13). 

In the reference orientation, the total intensity and cir- 
cular polarization parts of the beam response can be ex- 
panded in scalar spherical harmonics, e.g. 

i{e) =J2blim)nim){e), (22) 

Im 

where the sum is over I > and \m\ < I. The multi- 
poles bjj^-^ for the circular polarization are defined anal- 
ogously. For the beam response tensor Bat we must ex- 
pand in the transverse, trace- free tensor harmonics. Here 
we follow the coordinate-dependent approach of Ref. [13] 
(although for some applications the coordinate-free ap- 
proach of Ref. [32] is more convenient): 

B„6(e) = ^6f,„)F(L)„6, (23) 

Plm 

where the sum is over I > 2, \m\ < I, and the two types 
of transverse trace-free harmonics P = G (for Gradient, 
often called electric) or C (Curl, often called magnetic). 
All multipoles satisfy = (— l)'"6(;_m) since the fields 
are real and we have adopted the Condon-Shortley phase 
for the spherical harmonics. The Stokes parameters / 
and V for the sky can be similarly expanded in multipoles 
a^j^j and a^^,,^) , and the linear polarization in multipoles 

«(lm) and aC„^). 

The tensor harmonics are derived from the scalar har- 
monics by covariant differentiation over the sphere [13] 
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(see also the Appendix). Performing the differentiation 
gives [33] 



rGab 



1 



' (Im) = {-2Y(lrr,)m m + 2Y(^irr,)m* m*) , (24) 



Cab 
(Im) 



1 



-^{-2Y(^im)miS'm- 2YQm)Tn* i^m*) , (25) 



where m = [ae + ia^)/\/2. The ±2^(im) the spin- 
weight 2 harmonics: 



Ni 

±2Y(lm) = -^{W(im) ± iX(;„)), 

where iV; = [2(Z - 2)!/(/ + 2)!]V2, and 
52 



(26) 



-CSc2 6'^y(;„), 



Yi 



(Im) 



(27) 

(28) 



Explicit expressions for Vl^(/„j) and X(^i„-i^^ with the deriva- 
tives eliminated are given in Ref. [34]. Note that our 
convention for the spin- weight functions follows Ref. [29] , 
which differs from Goldberg et al. [35] by the inclusion 
of the factor (—1)"* in the definition of the spherical har- 
monics. 



A. Extracting the beam multipoles 

With the detector in the reference orientation, the 

spherical polar bases fixed relative to the sky and detec- 
tor coincide. We can extract the beam multipoles from 
the effective Stokes parameters on this polar basis using 
the orthonormality of the scalar and tensor harmonics 
(e.g. Ref. [13]). For the linear polarization, we have 

b'Cim) ± = 7f / T ^^)±2^(L) d^^, (29) 

which is the inverse of the expansion of Q and U in spin- 
weight 2 harmonics: 



-^(g ± iU) = ^(6|l„) T ibfim))-^2Y^imy 

Im 



(30) 



(The sum is over I > 2 and jmj <l.) Note that with our 
conventions for the polarization basis vectors, Q±iC/ is a 
spin-weight ^2 quantity, which differs from some authors 
(notably Ref. [14]). 
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FIG. 1. The orientation of the detector is specified by the 
three Euler angles {0, 0, V'} which takes the {fi} frame (i = x, 
y, z), which is fixed relative to the sky, onto the {<t^} frame 
which is fixed relative to the detector. The nominal main 
beam of the detector is along cr^, and the co-polar direction 
along the main beam is a'y. 

Some care is needed in extracting the beam multipoles 
at the north and south poles since the Stokes parameters 
on the polar basis are ill-defined there. However, using 
Eqs. (19) and (20), it is simple to show that any well- 
defined polarization field must have (p dependence going 
like Q ±iU (X cxp {zt2i(f>) on the polar basis at the north 
pole (9 = 0). This is consistent with Eq. (30) since the 
spin-weight ±2 harmonics satisfy 



±2Y^lm)=Sm^2\l^^^e^''^ 



(31) 



at ^ = 0. With finitely sampled, simulated data for 

£co and Eaoss, the contribution to hji„,-^ from samples 
on (or very near to) the north pole can be treated by 
approximating ±2y(;m) with Eq. (31), and absorbing the 
exp(±2i(/)) into {Q =F iU) exp(±2i0) which is then well- 
defined by the data. Similar problems occur at the south 
pole, but since the beam has virtually no power there 
this is not problematic. 



IV. BEAM ROTATION AND CONVOLUTION 

The kinematics of the experiment can be specified by 
a scan strategy which describes the rotation necessary to 
take each detector from the reference orientation to its 
orientation at the specified time in the scan. For simplic- 
ity we consider only a single detector, but our approach 
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could easily be generalised to experiments with multi- 
ple detectors. The rotation is specified by its Euler an- 
gles {(j), 6, V'}, such that first we rotate in a right-handed 
sense by tp about cr^, then by 9 about cry, and finally by 
(p about (7z again. We denote the rotation by D{<p, 6, tp) 
so that the image of CTj is D{^, 6, ip)cri = a'^ for i = x, y, 
z. For tp = 0, 

<j'^ = sin ^(cos (}f) <Ta; -I- sin <Ty ) + cos 9 az 

(Ty = — sirKpCTx + COS(p(Ty. (32) 

For non-zero ^p, cr'^ and a'y are additionally rotated by -0 
about a'z (see Fig. 1). 

Under the rotation D{(p, 9, tp) the beam Stokes param- 
eters I and V rotate as scalar fields so that e.g. 



i{e)^i[D{-i,,-9,-^)e], 



(33) 



where D{—tp, —9, —(p) is the inverse rotation. The beam 
response tensor Bab{e) rotates as a rank-two tensor, so 
that 



Bab{e) ^ Aa'^Ab^'Bcc, [D{-^P, -9, • 



(34) 



where Aa'^ is the SO (3) rotation matrix representing 
D{(p,9,ip) (see Appendix). Since we are describing the 
beams in multipole space, we must consider the transfor- 
mation properties of the scalar and tensor harmonics un- 
der the rotations given in Eqs. (33) and (34) respectively. 
The scalar spherical harmonics transform irreducibly un- 
der rotations as [36] 

Y{im){e)^ Dl,^{cP,9,^)Y^im'){e), (35) 

\m'\<l 

and, as we show in the Appendix, the same is true of the 
tensor harmonics: 

Y{im)ab{e)-^ ^ £'^,„(0,6»,V')y(im')ab(e). (36) 

\m'\<l 

Here, the Dl^,^{(p,9,'ip) are Wigner's £)-matrices. With 
our conventions for the Euler angles, we have 



D'^,^{,P,9,^P) 



-im' 4> il 



-imip 



(37) 



where 



mnK ) ' {l + m- t)\{l -n- t)\{t + n- m)W. 



X [cos(6'/2)] 



2l+m-n-2t 



sm{9/2)] 



n—m-\-2t 



• (38) 



The sum is over integers t such that the arguments of the 
factorials are non-negative. 

Performing the integral over the sphere in Eq. (13) 
is now straightforward in multipole space using the or- 
thonormality of the harmonics. The final result for the 
total power as a function of orientation of the detector is 



P 

(Im') 



Dl^,{cp,9,i,). (39) 



The sum is over / > with |m| and |m'| < I, where we 
defined the linear polarization multipoles to be zero for 
I = and 1. Our result for the dependence of the total 
power on orientation is quite general; we have made no 
assumptions about the beam profile and level of cross- 
polar contamination. Equation (39) is one of the main 
results of this paper. Note that the function Wtot (0, 9, tp) 
is expressed as a linear combination of the Z)-matrices, 
which form a complete set for expanding single-valued 
(square-integrable) functions on the three-sphere. 



A. Fast convolution algorithms 

The right-hand side of Eq. (39) can be evaluated 
rapidly by making only minor modifications to the al- 
gorithm developed recently by Wandelt and Gorski [27] 
for the case of an unpolarized detector. The key to the 
algorithm is to factor the rotation D{(p, 9, tp) as follows: 



D{cP,9,tP)=Di 



so we may write 



7r/2, -7r/2, 9)D{0, 7r/2, tP + 7r/2), 

(40) 



D^„,((/>,e,V)= E [i^™M('/'- 7^/2, -77/2,^) 

\M\<1 

xDW(0,7r/2,V + V2)]. (41) 

The advantage of factoring the rotation in this way is 
that now the Euler angles only occur in complex expo- 
nentials, and we only need evaluate <iJ,„„/(^) a,t 9 — tt/2 
[since d';^^^,{-9) = dl^'mi9)]- The full three-sphere of ro- 
tations can now be calculated with a three-dimensional 
Fast Fourier Transform. The dl^^,{9) can be computed 
with the accurate recursive method in Ref. [37] (which 
can be further enhanced by making use of the symme- 
tries of the d^^^,). To perform the convolution to a res- 
olution corresponding to multipoles Zmax, for all possible 
orientations of the system, requires 0(^4iax) operations. 
This should be compared with the 0(1^^^^) operations re- 
quired for a brute force computation in pixel space. For 
experiments such as the Planck mission, where Zmax is of 
the order of a few thousand, the saving is considerable. 
For many experiments, the approximate azimuthal sym- 
metry of the beam limits the sum over m' in Eq. (39) to 
Iw'l <C Zmax- Since Wtot is then a slowly varying function 
of tp, it is possible to sample the tp variation much more 
sparsely than for 9 and (p, which effectively reduces the 
operations count to 0{l^^) [27]. 
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V. AXISYMMETRIC BEAMS 

In this section we consider the hmiting form of the 
general result, Eq. (39), when the beam is approximated 
as being pure co-polarized (i.e. having no cross-polar 
contamination), and axisymmetric. Typically, these ap- 
proximations hold well across the main beam for a well- 
defined linear polarimetry system. However, the approx- 
imations are not valid for describing the response of the 
system in the far side-lobes, where the complete expres- 
sion, Eq. (39), should be used. To quantify the errors 
introduced in a given experiment by assuming an ax- 
isymmetric beam requires detailed simulation with the 
apparatus developed in the earlier sections of this pa- 
per. Techniques for propagating these errors to the band 
power estimates of the power spectrum in a total power 
experiment have been developed recently [26], but the ex- 
tension to polarized experiments must await the detailed 
development of the full analysis pipeline for polarized 
data. 

For a pure co-polar beam, Eaoss vanishes in Eqs. (18)- 
(21). If we further assume axisymmetry, then £co is a 
function of 6' alone. Writing |£coP = B{9'), we have 
i = B{9'),V = 0, and 



Q±iU = -B{e')e 



±2i<j,' 



(42) 



It follows that b^i^) = unless |m| = 2. Furthermore, 
using Eq. (29) with m = 2, we see that 6^2) = *^(i2) 

hence ^(j_2) ^ ^''^[i-2)' ^^'^ beam response tensor for 
linear polarization, Babi is now fully specified by h'^i2y 
If we use Eqs. (27) and (28) to express the spin-weight 

harmonics in terms of Legendre functions, we find the 

G . 
(Z2)- 



following expression for hf 



irNf 121 + 1 



V2 



+ 



2{l-l)x--l{l-l){l-x'')-{l-4) 



(43) 



where primes denote differentiation with respect to a; = 
cos 9', and Pi {x) is the Ith Legendre polynomial. Note 
that 6|]2) is real, so 6^2) is imaginary. For the intensity 
multipoles we have = unless m = 0, with 



2tt 



21 + 1 
in 



Pi{x)B{x) dx. 



(44) 



Given the restrictions on the beam multipoles for an 
axisymmetric, co-polar beam, it is possible to simplify 
Eq. (39) for the total power received. Using the rela- 
tion [35] 



47r 



D^^_^{<j>,9,i;) = {-iy^-^,Y^^^^{9,<i>)e^^^ (45) 



between the _D-matrices and the spin- weight s harmonics 
for integer I, m, and s, and the reality of 6^^q^ and ^(^2)' 
Eq. (39) can be written in the form 

Wtot (e, (/>, V) « ^ [les - Qes cos 2V' + f/eff sin 2^ ]. (46) 
Here, 7eff is the usual beam smoothed intensity: 

/eff(0, </.) = 5] Wiali^^Y^i^){9, <P), (47) 



where the window function 



(48) 



Similarly, Qeff and Uea are beam smoothed Stokes pa- 
rameters on the spherical-polar basis, given by 

Im 

where, following Ref. [29], we have introduced the spin- 
weight 2 window function 



2W1 



-2V2 



2; + r('2)- 



(50) 



Note that Eq. (46) shows that the signal obtained by con- 
volving the pure co-polar, axisymmetric beam with the 
sky is equivalent to the response of an idealised co-polar 
detector, with vanishing beam width, on the smoothed 
sky. This result does not depend on any assumptions 
about the angular size of the beam response; for the po- 
larized contribution it is a consequence of the definition 
of the co-polar vector field, Eq. (16), which is the obvious 
generalisation of a constant vector field to the surface of 
sphere. Equations (46) (50) provide a complete descrip- 
tion of the power received in polarimetry experiments 
with axisymmetric, co-polar beams. 



A. Gaussian beams 

It is often the case that the axisymmetric beam profile 
is approximately Gaussian: 



5(6*') = -Bexp[-(1 - cos6'' )/(t2]. 



(51) 



where cr is a measure of the beam width. Note that we 
follow Bond and Efstathiou [38] in taking the beam to 
be Gaussian in 2sin(0'/2) rather than 9' . The former 
allows us to derive simple analytic results valid for any 
a. However, in most cases a Gaussian profile is only 
appropriate close to the beam axis, in which case the two 
definitions are almost indistinguishable. For the beam 
profile in Eq. (51), the window functions Wi and 2W1 
can be evaluated analytically: 
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Wi = 47rBe-"i;(a), (52) 

and 

2W1 = 27rB7Vf e-«{2(-l)'[(Z + 2){l - 1) + 6a]e-« 

+ [{f - 4.a)il - If + 12a^]ii{a) 

+ 4a(Z2 + Z + l-3a)i(_i(Q)}, (53) 

where ii{x) is a modified spherical Besscl function, and 
a = l/cr^. For cr^ <C 1, which is always the case for 
high resolution experiments, the window functions are 
very well approximated by their asymptotic expansions, 
which give 

Wi!=i2TrBa^cxj>[-l{l + l)a^/2], (54) 
2W1 « 2wBa'^ exp{-[l{l + 1) - 4](T^ /2}, (55) 

in full agreement with Ng and Liu [29]. For Z 3> 1, the 
polarized and unpolarized window functions in the small- 
scale limit are approximately equal, which can easily be 
verified by making a flat sky approximation. 

VI. CONCLUSION 

In this paper we have presented a multipole method 
for describing the response of an arbitrary detector and 
feed system. Wc fully include the effects of non-zero 
beam size, asymmetric beam patterns, and cross-polar 
contamination. Such inclusions are essential for the accu- 
rate modelling and interpretation of precision polarime- 
try data, such as that expected from the Planck mission. 
Working in multipole space, we derived a simple expres- 
sion for the response of the system, when convolved with 
the sky, as a function of the three Euler angles needed to 
describe a general orientation of the system. Given the 
form of this expression, it is straightforward to modify 
the fast algorithm of Wandelt and Gorski [27] to com- 
pute the system response for the entire three-sphere of 
orientations in 0{l^^) operations. Finally, we showed 
how, for the case of a pure co-polar, axisymmetric beam, 
the response can be described by spin-weighted window 
functions. This extended the results of Ref. [29] to ar- 
bitrary size beams, and we gave the exact form of the 
window functions for a beam with a Gaussian profile. 
Although our discussion has been in the context of CMB 
polarimetry experiments, the formalism introduced here 
should be useful in any applications involving anisotropic 
filtering of tensor fields on the sphere. 

The techniques described in this paper have now been 
implemented in the simulation pipeline for polarized 
channels of the Planck mission. 
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APPENDIX A: ROTATING THE TENSOR 
HARMONICS 

In this appendix we establish the transformation prop- 
erties of the transverse, trace- free tensor harmonics under 

active rotations. The action of the rotation D{(j>,9,tp) 
can be represented by a rotation matrix A";, such that 
an arbitrary vector u° rotates to v'" = A"bV^. Orthogo- 
nality of the rotation implies A'^i,Aa'^ = so the inverse 
rotation is w° = v'^Ab"-. Under an active rotation of the 
beam, the beam response tensor rotates to B'^^{e'') which 
is obtained by forward rotating the original tensor eval- 
uated at the back rotated position: 

BUel = A/^A,=^S,,,,(e'^A/). (Al) 

Under this transformation, the transverse nature of the 
tensor field ^06(6*^) is preserved. We shall demand the 
same transformation properties for the tensor harmonics 

yP 
{lm)ah' 

It is convenient to write the tensor harmonics in terms 
of covariant derivatives on the sphere of the scalar har- 
monics (e.g. Ref. [13]): 

^(L)a6 = (VaV,Y(,„,) - \9alN^Y(lm^ (A2) 
y^m)ab = Y(^"bV„VcF(;„) + ~e\V ,V .Y^i^)) , (A3) 

where Va is the covariant derivative on the unit sphere, 
Qab = gab ~ ^a&b IS the (iuduccd) mctric (with gab the 
metric of Euclidean 3-space) and tab = ^abcS-'^ is the pro- 
jected alternating tensor. The covariant derivative Va is 
obtained from the 3-dimensional (flat) covariant deriva- 
tive Va by total projection: 

Van...c=~9i~9l---~9l^dTe...f, (A4) 
for an arbitrary tensor Ti,,„c- Making use of the results 

(Val^(';H)le= = K\^bY(lm))UK,^ (A5) 
(VaV6y(';„))|ec = Aa'^^A(,''HV<i,V<i,F(i„))|e.3A,3c, (A6) 

it is straightforward to prove that the rotated tensor har- 
monics arc obtained by replacing the scalar harmonics 
by their rotated counterparts Yli^^ie-'^) = Y(^im){e'^ Aif) 
in Eqs. (A2) and (A3). Since the Zth scalar harmonics 
transform irreducibly under rotations as (e.g. Ref. [36]) 

Ylim){^) = E Dl,^{'PM)Y(im'){e), (A7) 

\m'\<l 
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the tensor harmonics inherit the same transformation 
law: 

^(L) (e) = E ^'m'm{4>, 0, miL')ie)- (A8) 

\m'\<l 
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